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ABSTRACT 

Aims. The aim is to study structure and polarization properties of the stars and planet systems in the active Orion H II-region. 
Methods. We performed AO-assisted high-resolution imaging polarimetry on selected Orion proplyds close to the Trapezium stars in 
the J, H, and K bands. Differential polarimetric images of one of the largest and brightest proplyds are interpreted using 3D radiation 
transfer simulations based on the Monte Carlo method. 

Results. Although not fully resolvable by ground-based observations, the circumstellar material can be mapped with polarimetry. 
We present constraints on the disk parameters of the giant proplyd 177-341. We tested whether dust models with different grain size 
distributions could explain the observed extent of the polarization patterns and find that simple models with larger grains will not 
reproduce the spectral polarization behavior. 

Conclusions. The technique of polarimetric differential imaging (PDI) in the NIR provides a good opportunity to study the structure 
of the Orion proplyds. 

Key words. Polarization - ISM: HII regions - Stars: 177-341 - Stars: early-type - Stars: circumstellar matter - Stars: planetary 
systems: proto planetary disks 



1. Introduction 

The Orion molecular clouds at a distance of about 460 par- 
sees are the nearest location of ongoing massive star formation. 
This region also harbors one of the densest young stellar clus- 
ters dMcCaughrean & StauffeJll994l) . the Orion Nebula Cluster 
(ONC). The core of the ONC is known as the Trapezium clus- 
ter, named after the four massive OB stars at its very center. The 
dominant Trapezium star is 8 l Ori C, which generates most of 
the intense ultraviolet radiation field that photo-dissociates and 
photo-ionizes the surrounding material. Th e discovery of proto- 
planetary disks (proplyds) in this region bv lO'Dell etai1(ll993l) 
with the Hubble Space Telescope (HST) has provided an un- 
precedented opportunity to study planet-forming disks around 
young stars. These objects are surrounded by bright, crescent, 
and tear-drop shaped ionization fronts (IF) resulting from gas 
flow off of the disk. At greater distances and shielded from radi- 
ation, several proplyds without ionization front have been ob- 
served. They are seen in silhouette against the nebular back- 
ground dBallv et al.ll2000l) . 

New insights into the nature of proplyds have been ac- 
quired through hig h-resolution multi-spectral imaging with HST 
(Ball y et all [l 9981) . follow -up imaging, and spectroscopy at a 
variety of wavele ngths from the ground, and intensive theo- 
retica l modeling dJohnstone et alj 119981: IStorzer & Hollen bach 
1999, e.g.). As a result, an irradiated disk model has been de- 
veloped to account for the presently known observed features. 

Near infrared (NIR) observations of the Trapezium clus- 
ter dHillenbrand et all 1 1 9981: lLada et al.1 [20041) have shown that 
most young stars have excess emission in this wavelength 
regime. This is interpreted as the presence of circumstellar disks. 
According to the study bv lLada et al.l (|2000), 80% of the stars in 



the Trapezium cluster show an infrared excess and are likely to 
be surrounded by disks. 

Currently much activity is devoted to obtain physical prop- 
erties of the circu mstellar dust distribution and these disks. 
lOlczak et al.l d2006) showed that encounters of stars and circum- 
stellar disks in the Trapezium cluster are quite frequent and thus 
can cause a considerable contribution to the mass loss and the 
truncation of proplyds. Encounter-triggered planet formation is 
also possible a s an impo rtant supplement to triggering by super 
massive stars (Ballv 2003) and supernovae-triggering of the OB- 
stars (Hester & Desch 2005). The disks of the proplyds close to 
the Trapezium with a typical radius of 20-80 AU are not resolv- 
able directly. In the Orion nebula the YSO 177-341, also known 
as HST1 is one of th e best-studied proplyds because of its large 
size and orien tation dBallv et al.lll998b iHennev & Q'Delllll999l: 
lHenney|2000l) . The very bright ionization front outshines details 
of its inner structure. 

iMonin et al.l d2006) showed that polarimetry can be used 
to trace the relative orientation of disks in young binary sys- 
tems. The polarization structure of a massive young stellar object 
(YSO) in OMC-1 reveal ed a surrounding disk/bipolar outflow 
system (Jiang et al. 2005). Near-infrared polarimetry is a proven 
remed y for mapping hidden effects and structur es (see, for ex- 
ample, Ageorg es et al.|[l997l iMever et al.l 12006b . In the case of 
dust around a bright star, only the scattered light from the disk 
is expected to be polarized, the unpolarized light from the cen- 
tral star is suppressed by imaging polarimetry. Therefore this 
technique can improve the contrast and the detection-threshold, 
and it provides a possibility to map struct ures very close to 
bright stars wi thout coronographic masks dKuhn et alj 12001; 
lHales et al.l2006l) . Polarimetric images can also trace the sources 
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Fig. 1. Near-Infrared (left, KHJ = RGB) image of the giant pro- 
plyd 177-341 and its neighbors in comparison to the o ptical HST 
image (right, red is [Nil], green is H ff , blue is [Q Il. lBallv et al.l 
1998), with north facing upwards and east to the left. The pro- 
plyd 178-342 is covered by the ionization front in the optical, 
but barely visible at the end of the tail. 



of illumination of the r eflection nebula in star-forming regions 
dHashimoto et alJl2007h . 

We carried out a JHK polarimetry survey of selected pro- 
plyds in the Trapezium Cluster to study dust and disk parameters 
in this dense star-forming region. Those disks of the proplyds 
close to the Trapezium that have a bright IF are so far only seen 
as IR excesses. With the polarimetric images we want to map 
the circumstellar dust distribution and extract parameter ranges 
for the disks to classify the Trapezium proplyds. In this paper we 
want to present the results from our first polarimetric study of the 
Orion proplyds. We begin in Sect. 2 with the description of the 
observational data and summarize in Sect. 3 the basic properties 
of one of the sources extracted from the polarized data set. In 
the next section we consider different analytical techniques for a 
more detailed view of the proplyds properties including model- 
ing of circumstellar dust particles and qualitative comparison to 
simulated NIR polarimetric images. A summary and an outlook 
of forthcoming surveys is noted in Sect. 5. 



2. Observations and data reduction 

High-resolution near-infrared polarimetric images were ac- 
quired in 2001 with the ESO Adaptive Optics System ADONIS 
and the SHARP 11+ cam era at the La Silla 3.6v$\ telescope 
dAgeorges & Walshlll999h . A rotating wire grid polarizer pro- 
vided polarization angles from 0° to 150° with a 30° incre- 
ment. The sources are 10 proplyds close to the Trapezium stars, 
embedded in three pointings with a field-of-view (FOV) of 
8.96"x8.96" and a 0.035" plate scale. Images were taken over 
three nights in the J, H, and K band. No occulting masks were 
used in these observations due to multiple sources per FOV. The 
chosen detector integration time for a good signal-to-noise ratio 
(SNR) and dynamic range was 15 s per polarization angle and 
FOV. The total time spent on each source was 6750 s, with sev- 
eral measurements per angle in different nights (375 s per source, 
angle, and band). 9 l Ori C was used for wavefront sensing of the 
adaptive optics. The natural seeing varied from 0.5" to 1". The 
achieved Strehl ratio in JHK was 6%, 13%, and 20% for the most 



distant FOV. Off-source sky for each angle and dark frames were 
taken, as well as frames with unpolarized reference stars. 

The data reduction was carried out using the 
DPUSER software for astronomical image analysis 
(T Ott, ht tp://www.mpe.mpg.de/~ott/dpuser/, see also 
lEckart & Duhouxf 1 19901) and IDL routines. Due to some 
cooling problems of the detector, some part of the data had 
to be sorted out by a statistical filter algorithm followed by a 
manual inspection. The raw images were sky-subtracted and flat 
field-corrected. 

The non-polarized light consists mainly of the central star's 
light and the uncorrected AO residuals, which is the domi- 
nant noise source of ground-based AO-observations in the NIR. 
Polarimetric images are very sensitive to minor deviations of 
shift-and-add operations. Therefore this reduction step was care- 
fully carried out and manually optimized for every single frame. 
A minimal number of three frames (spanning more than 90° po- 
sition angle) is required to determine the linear polarization and 
its position angles. The redundant data set of three orthogonal 
polarization component pairs increases the SNR and reduces the 
instrumental polarization. By fitting the oversampled polariza- 
tion curve with a cosine function, the linear polarization angle 
and degree maps (YIl) have been determined. Polarimetric dif- 
ferential images (PDI) can be obtained using Qo = 7o° - ho- 
(Stokes Q-parameter), Q 30 = ho° - hw, and £>60 = ho° - hsw- 
In the widely used notation of linear polarized light this converts 
to Q - Qo, and 

Q<p-Q- cos(20) 



U = 



with d> = 30° or d> = 60°, and 



sin(2</>) 



n, = 



Ve 2 + u 2 



1 u 

= - arctan — 

2 Q 



ESO program-ID: 066.C-0219 



The Q- and U-parameters are used for better comparison to the 
simulated images of the model. 

3. Results 

Some of the sources show a clearly extended structure in the 
NIR referring to an illuminated dust envelope. The surround- 
ing nebula emission is not as strong and spatially varying as in 
the optical. The properties of the targets are listed in Table Q] 
including the IF-sizes and NIR magnitudes, J-H and H-K col- 
ors. These were photometrically calibrated with the standard star 
HD98161 from the 2MASS catalog. This unpolarized star was 
also used for polarimetric calibration. The aperture polarization 
of the targets is shown in Table [2] Also the integrated polariza- 
tion in two annuli is shown for comparison and classification 
of circumstellar material. The inner regi on is usually domi nated 
by alignment and subtraction residuals dHales et al.l l2006). The 
annulus of 0.3"-0.5" refers to the circumstellar dust, which is il- 
luminated by the central star. The outer annulus should be dom- 
inated by polarized light from the surrounding reflecting nebula 
where the source of illumination is unknown. 

The polarization values as a function of the wavelength are 
plotted in Fig. [2] Although the integrated polarization degree and 
angle is highly affected by minor deviations in shift-and-add op- 
erations and PSF effects (see Sect. 14.21 ), the average values give 
a good overview of the circumstellar polarization distribution. 
The nebula emission in the outer annuli is polarized from 10% 
in K-band to 20% in J-band. This is in good agreement with 
a recent wide-field polarization study of the ONC in the NIR 
dTamura et al.ll2006h . 
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Fig. 2. Spectral distribution of polarization for J, H, and K bands. 
In an aperture of 0.1" the light of the central star the scattering 
by foreground material dominates. In an annulus close to the star 
the polarization is higher due to the circumstellar envelope and 
disk. The outer annulus is dominated by the reflection nebula and 
is highly polarized. The dashed line is the average polarization 
spectrum, the dotted line is 177-341. 
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Fig. 3. Polarimetric differential images (KHJ = RGB) of the 
sources in the FOV shown in Fig. \T\ Top: Proplyd 177-341 in 
Qo (left) and £>60 (right). A clear butterfly pattern is visible in 
all bands. Mid and Bottom: 178-342 and 178-343 show a tre- 
foil that is supposed to be an AO-effect. With this comparison of 
sources in one FOV, lens- (e.g. astigmatism) and PSF-effects can 
be excluded as the origin for the results of 177-341. 



To remove all non-polarized light from the sources, we sub- 
tracted orthogonal components from each other. These images 
of the Stokes Q- or U-parameters show a characteristic trefoil 
pattern for most of the sources (Fig. 0). The trefoil is a mode 
in the AO wavefront reconstruction that may contain residual 



power if the AO-loop is run under variable atmospheric condi- 
tions. However, for the giant proplyd 177-341 we find significant 
butterfly patterns or two lobes, which is the polarization pattern 
expected for a circumstellar dust envelope. 

The sources around 177-341 (Fig. [3]) lie in a relatively dark 
region of the nebula. This is why we concentrate here in particu- 
lar on the field of the proplyds around 177-341. The NIR-image 
in Figure Q] shows quite a different vie w of the proplyds mor- 
phology than the known optical images (Bal ly et al.lfl998l) . First 
of all the bright ionization front, which covers inner structures 
of the proplyds in the optical is not visible in the NIR-images. 
This reveals, for example, a previously unknown source behind 
the tail of the proplyd 177-341 (Fig.Q] left). If the new source 
is a proplyd, it would be n amed 178-342, followi ng the designa- 
tion system introduced by|0 'Dell & Wed Jl 994). The observed 
JHK-colors differ from the neighbor 177-341, 178-342 seems 
redder (Table[TJ. This may be due to extinction by the overlying 
ionization front. 

In Fig. [4] the polarization vectors patterns of 177-341 are 
shown. Close to the central source in the dust envelope, a cen- 
trosymmetric pattern is visible in relatively low polarization in 
the range of 5-10%. The reflection nebula in the outer regions 
of the proplyd shows a uniformly oriented polarization pattern 
with the same direction in all bands. This is almost perpendicu- 
lar to the direction of 6 l Ori C, which is most likely the source 
of illumination of this part of the nebula. 



4. Model and discussion 

4.1. Polarization properties 

Circumstellar accretion disks mark optically thick regions of the 
reflection nebula where the polarization patterns clearly depart 
from centrosymmetry. The vectors in these regions usually show 
low ampli tudes of polarization of about 10% due to multiple 
scattering (Bas tien & Menardl 19881) . In the aperture of 0.35" for 
J-band to 0.5" for K-band, the centrosymmetric pattern is visible 
to the left in Fig. [4] A region of shortened and aligned vectors is 
not clearly resolvable in this view. 

The outer regions show higher polarization and an aligned 
pattern from the illuminated reflection nebula from 9 1 Ori C. The 
extent of the polarization pattern around 177-341 is determined 
in the radial plots to the right of Fig. [4]as follows. The fraction of 
polarized light as a function of the distance from the star shows 
a plateau of low polarization of the circumstellar envelope and 
is increasing for higher radii due to the nebula polarization. For 
a perfect centrosymmetry, the profiles along the position angle 
of the Q- and U-images should resemble a cosine. The ampli- 
tude of the cosine fitted to the Q- and U-profiles as a function 
of the distance then clearly defines the extent. The maximum ra- 
dius increases from J- to K-band and a small plateau around the 
maximum extent is visible in all bands. 

For a closer look at the spatial proportions, we examined the 
profile of the Q- and U-images by adding 5-degree bins over 
the determined radial range. These azimuthal Q-profiles show 
the characteristic cosine or W-curve for a dust envelope (Fig. [4] 
middle). In our data the profiles are not symmetric around zero as 
they should be, and this is equivalent to an offset in polarized flux 
and is related to effects of the PSFs. The long distance of about 
26" from 177-341 to the AO guide star 1 Ori C (Table CD re - 
sults in elongated PSFs for the observations dCresci et alj |2005). 
The direction of elongation is towards 9 l Ori C (134° PA.). A 
shifting of the profiles can be seen, after convolving simulated 
polarization images with the PSF of the observations (Figs. [13] 
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Table 1. JHK-magnitudes and colors. 



Designation 


R.A. 


DECL. 


Dist. to 


IF a (opt.) 


Magnitude b 




Color 




(J2000) 


(J2000) 


6 l Ori C 


Tail 


J 


H 


K s 


J-H 


H-K 


157-323 


05 35 15.71 


-05 23 22.59 


11.25" 


<0.7" 


9.88 


9.05 


8.93 


0.83 


0.11 


158-323 


05 35 15.82 


-05 23 22.50 


9.75" 


1.8-2.5" 


9.73 


8.83 


8.53 


0.91 


0.3 


158-326 


05 35 15.83 


-05 23 25.62 


10.2" 


0.85" 


12.9 


11.34 


10.65 


1.56 


0.69 


158-327 


05 35 15.79 


-05 23 26.61 


10.53" 


1.0-1.5" 


11.68 


10.4 


9.66 


1.28 


0.75 


163-317 


05 35 16.27 


-05 23 16.72 


6.42" 


1.6" 


10.14 


9.41 


9.2 


0.72 


0.21 


166-316 


05 35 16.60 


-05 23 16.32 


7.36" 


0.4" 


11.4 


10.47 


10.2 


0.93 


0.27 


167-317 


05 35 16.73 


-05 23 16.63 


7.08" 


0.67"x2.1" 


10.37 


9.41 


8.94 


0.96 


0.47 


177-341 


05 35 17.67 


-05 23 40.96 


25.98" 


0.8"x3.5" 


12.44 


11.83 


11.53 


0.61 


0.3 


178-342 


05 35 17.76 


-05 23 42.50 


27.2" 




12.93 


11.67 


11.2 


1.26 


0.46 


178-344 


05 35 17.78 


-05 23 44.25 


29.1" 




12.083 


11.28 


11.1 


0.8 


0.17 


HD98161 


11 17 12.01 


-38 00 51.72 






6.075 


6.019 


5.992 


0.056 


0.027 



a lBallv etal.lfi99"8l) 

b magnitude errors ~ 0.1 



Table 2. Integrated polarization of the observed proplyds and surrounding reflection nebula in an aperture of 0.1", in an annulus of 
0.3"-0.5" and 0.7"-1.0", respectively. 





p.j 




Aperture 0. 1 " 






Annulus 0.3"-0.5" 






Annulus 0.7"- 1.0" 








Pol. 


t <°\ 






Pol.[%] 




■t j 


Pol.[%] 




1 3 1 

■L J 


157-323 


J 


18.11 


±0.70 


35.5 


±1.0 


14.38 


±0.25 


45.5 


±0.3 


14.89 


±0.37 


38.7 


±0.8 




H 


4.76 


+0.26 


84.4 


±5.0 


6.41 


±1.96 


22.7 


±0.2 


11.95 


±0.68 


26.0 


±1.0 




Ks 


3.71 


±1.91 


62.1 


±4.7 


2.77 


±0.64 


-67.2 


±29.6 


11.25 


±0.71 


-56.7 


±0.9 


158-323 


J 


19.61 


±2.09 


21.0 


±2.1 


13.33 


±1.65 


44.6 


±0.3 


12.11 


±0.77 


40.3 


±0.5 




H 


3.58 


±0.18 


61.3 


±2.7 


5.33 


±2.22 


20.8 


±1.8 


9.30 


±0.63 


23.3 


±0.4 




K.v 


4.12 


±1.22 


24.4 


±6.0 


2.72 


±0.92 


-60.4 


±31.4 


10.31 


±0.37 


-51.5 


±0.6 


158-326 


J 


45.32 


±4.52 


35.9 


±2.1 


32.53 


±0.71 


29.9 


±0.7 


46.87 


±1.53 


29.4 


±0.1 




H 


9.12 


±2.50 


89.2 


±1.6 


9.26 


±1.31 


29.9 


±0.6 


20.06 


±0.74 


29.1 


±0.1 




K, 


3.08 


±1.17 


-76.4 


±10.5 


2.61 


±0.46 


42.2 


±7.4 


4.30 


±0.61 


17.6 


±2.0 


158-327 


J 


19.14 


±3.42 


32.8 


±0.4 


17.27 


±0.93 


40.9 


±1.1 


18.24 


±0.85 


35.0 


±0.2 




H 


6.94 


±1.69 


-78.4 


±3.4 


7.95 


±1.25 


34.0 


±0.8 


11.69 


±0.21 


33.7 


±0.3 




K, 


3.07 


±1.84 


-40.2 


±16.0 


3.76 


±0.42 


58.1 


±7.4 


1.93 


±0.22 


-88.7 


±2.3 


163-317 


J 


24.77 


±2.72 


34.9 


±1.5 


9.92 


±0.44 


59.5 


±2.9 


13.23 


±1.65 


71.0 


±0.3 




H 


2.89 


±1.66 


-29.5 


±23.1 


6.41 


±1.91 


68.3 


±4.3 


12.57 


±0.87 


-87.3 


±0.3 




Ks 


3.96 


±1.61 


-79.1 


±5.5 


5.10 


±0.98 


29.9 


±1.0 


3.12 


±0.43 


-24.8 


±21.0 


166-316 


J 


25.50 


±1.58 


44.3 


±1.3 


14.28 


±1.95 


54.8 


±1.7 


9.75 


±1.33 


77.5 


±0.5 




H 


3.00 


±0.94 


36.2 


±8.2 


7.46 


±0.71 


69.5 


±2.8 


5.88 


±0.29 


76.4 


±1.1 




K.s 


5.35 


±2.07 


-56.6 


±33.3 


2.22 


±1.09 


82.4 


±9.8 


4.09 


±0.34 


-44.2 


±35.1 


167-317 


J 


23.94 


±1.38 


30.8 


±1.7 


11.70 


±1.25 


38.5 


±1.0 


15.28 


±2.48 


11.9 


±2.3 




H 


5.86 


±1.48 


-10.8 


±1.5 


7.84 


±2.30 


56.8 


±3.0 


13.34 


±6.57 


21.4 


±6.4 




K, 


1.14 


±0.75 


85.6 


±11.5 


5.09 


±0.45 


68.4 


±2.9 


8.26 


±6.54 


19.1 


±9.6 


177-341 


J 


5.54 


±1.72 


-17.0 


±4.5 


9.68 


±1.64 


-6.1 


±3.8 


20.38 


±2.59 


-3.6 


±0.9 




H 


1.64 


±0.96 


36.5 


±4.1 


10.27 


±0.80 


27.3 


±1.7 


19.90 


±2.02 


17.4 


±0.5 




K, 


3.42 


±0.20 


30.3 


±0.3 


0.80 


±1.06 


-46.3 


±33.0 


11.14 


±0.70 


-58.0 


±0.4 


178-342 


J 


17.12 


±3.77 


50.7 


±0.6 


20.48 


±2.03 


15.4 


±1.6 


25.18 


±2.87 


7.5 


±0.2 




H 


6.30 


±2.07 


-41.3 


±35.5 


11.02 


±3.93 


16.2 


±4.8 


23.94 


±1.99 


6.4 


±0.6 




Ks 


11.11 


±3.19 


63.5 


±4.4 


0.36 


±1.14 


87.5 


±30.6 


3.14 


±1.08 


-48.2 


±34.2 


178-343 


J 


11.08 


±0.91 


59.2 


±4.3 


17.64 


±2.59 


25.8 


±1.6 


23.57 


±1.01 


15.4 


±0.5 




H 


11.94 


±3.99 


-46.3 


±0.4 


11.94 


±3.03 


23.5 


±4.2 


24.00 


±3.99 


8.3 


±0.6 




K s 


10.14 


±3.25 


87.2 


±0.5 


2.35 


±0.73 


-39.6 


±25.6 


8.27 


±1.35 


-62.6 


±0.6 



and [14] only available electronically); therefore, the profiles of 
the convolved images of the data and simulations were shifted to 
be symmetric around zero so as to achieve a result closer to the 
actual centrosymmetric polarization pattern. This can be seen as 
a bootstrap calibration of the polarimetric images on the basis 
of the centrosymmetry of the envelope polarization. This has an 
impact on the pattern of the surrounding reflection nebula, which 
is highly polarized, but the total flux is low. The relative changes 
in the region of centrosymmetry are irrelevant due to the higher 
flux, while the analysis of the individual Q- and U-images in this 
study is not affected. 



Since the Q- and U-images are not acquired simultaneously, 
there are deviations due to residual speckles and changes in see- 
ing in addition to the elongated PSFs. These errors are empha- 
sized during the determination of fl^ = VG 2 + U 2 . As a con- 
sequence, no clear features are seen in the polarization degree 
images, and no direct comparison to convolved simulation im- 
ages is possible. Hence, we tried to deconvolve the observational 
data in this case. For the PSF estimation the F OV with 177- 
341 w as analyzed with the StarFinder program (Diolait ret al.l 
2000). Three point-like sources around 177-341 were used to 
extract a field PSF (Fig. Q3] only available electronically). For 
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Fig. 4. The polarization vector maps of 177-341 show a significant pattern. Left: Stokes-Q-images with intensity contours (logarith- 
mic, dotted line is 50%) in J-, H-, and K-bands (from top to bottom), the vector length is proportional to the degree of polarization. 
The characteristic butterfly pattern refers to an extended polarized source around the central star. The reflection nebula shows a 
uniformly oriented polarization vector pattern, which is almost perpendicular to the direction of 6 l C Ori (indicated by the arrow 
in the lower left corner). For comparison, the PSF (with dotted 50% contour line) is plotted on the same scale in the lower right 
corner. Right: Radial polarization profile, the linear polarization degree increases with the distance to the central source. The am- 
plitude of the cosine-fit of the Q- and U-profiles (along the RA.) over the radius is also plotted, which clearly defines the extent of 
the polarization pattern. The ranges are plotted with vertical lines. Middle: Stokes-Q-profile against the angle in Jansky: 5 degree 
bins in the Q-component against the rotation angle in the range determined with the radial plot to the right and shown in the two 
Q/U-doughnut-plots (0°-angle facing to the west, rotation counterclockwise). A slightly greater amplitude in the range of about 
50°-230° (interval II) of the Q-profile is visible in all bands. The ratio of the fitted amplitudes in II and 12 (Al, A2) are also plotted 
as a function of the radius to the right in linear range from to 2. 



the deconvolution of the i ndividual in tensity images 1$, the it- 
erative Lucy-Richardson dLucvl IT974) algorithm was chosen, 
which is known to be robust and stable even under difficult 
SNR conditions dPruksch & Fleischmann| [l998'). For the beam 
reconstruction, the deconvolved /^-images were convolved with 
a Gaussian of 0.25" FWHM, before the polarization degree im- 
ages were generated. The diffraction limit at K-band is 0.15" for 
the 3.6m telescope, but under moderate seeing conditions and 
these low Strehl ratios of the observations (Sect. |2|, an empir- 
ical li mit for the achieved r esolution is 1.5-2 diffraction limit 
units dBonaccini et alJ [T997). Without this smoothing, the de- 
convolved images are also dominated by high-frequency noise 
features enhanced p rogressively with the number of iterations 
dLanteri et alj Il999f) . The algorithm converges fastest at the 
brightest source structures dLucvll 19741) . Therefore choosing too 



few iterations results in images of undefined resolution and qual- 
ity. Fr om our extensive experience with the algorithm dOtt et al.l 
1999), we chose 10 3 iterations for the case of a single star dom- 
inating a single circumstellar structure. 

For the position angle of the disk, we compared the HST 
optical images and simulations to the polarization degree maps, 
obtained from deconvolved intensity images (Fig. [5). The pa- 
rameters of the simulation are discussed in Sect. 14.21 but here 
a model with highly polarized envelope is shown to illustrate 
the disk signature in the envelope signal. In the optical images 
a small portion of the disk is seen in silhouette around the cen- 
tral star. Almost the same RA. is found as an elongated region 
of low polarization in the NIR polarization images and can be 
quantified to ~ 144°. The diameter of the disk is about 300AU 
in this view. It must be stated that the direction of elongation of 
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Fig. 5. Left: HST optical images dBallv et alj2000l) of 177-341 through the filters f656N (H a, top) and f502N ([O III], bottom). The 
images have been subtracted by a Gaussian-smoothed version of itself, which enhances the local contrast and attenuates extended 
emission. A possible disk orientation is marked. Right top: Polarization degree images of 177-341 in J-, H-, and K-bands (from 
left to right). The individual intensity images 1$ were deconvolved before the calculation of the polarization degree, which is very 
sensitive to PSF-efFects. Right bottom: Polarization degree images for a MC simulation with disk inclination of 75°. 



the PSFs differs only 10° from the obtained PA. and the elon- 
gated low polarization region could just be an artefact. But the 
extent of the region of 300AU clearly exceeds the PSF dimen- 
sions at least in the K-band. The PA. is also confirmed by the 
following analysis. 

There is a noticeable modulation of the amplitudes of the 
profiles, which can constrain disk parameters if we consider the 
optically thick accretion disks. In regions where the disk lies in 
front of the illuminating source from the observer's point of view 
the amplitude should be attenuated. The same applies where the 
disk blocks the propagation of the central light to outer regions of 
the envelope. This can be associated with the amplitude modu- 
lation (AM) technique used for communication and signal trans- 
mission. The envelope scattering is the carrier signal, i. e. a co- 
sine with fixed phase in the Q- and U-profiles. Since there is 
no reference amplitude and we have only two cycle periods (see 
Fig. 0J, only the ratio between the two can be determined as 
the modulation signal. To investigate these effects, we generated 
amplitude ratio maps, which are defined in polar coordinates as 
follows. For every position angle and radius (ro, <po) the Q-profile 
over a range from (r ,0o-9O°) to (r, 0o+9O°) is fitted with a func- 
tion f(<p) = Aq + A cos {2(f> - (/)q) with free amplitude A and am- 
plitude offset An, but fixed phase offset 0q = 0° and <pu — 45°, 
respectively. The best-fitting amplitude is divided by the one of 
the remaining 180° (for clarification see Fig. [4] middle: dash- 
dotted lines). The cosine-fitting incorporates the model of the 
fixed carrier signal amplitude. The individual fitting over the ra- 
dius r allows arbitrary radial density and polarization profiles for 
the envelope, but rotational asymmetries like flattened envelopes 
are not considered. The result has an advantage over the pure Q- 
/U-profiles, since the attenuation of the signal is determined over 
the whole position angle, and the SNR is increased under the as- 
sumption of centrosymmetric envelope profiles. In principle the 



achievable spatial resolution is also enhanced by the high sen- 
sitivity of the signal to modulations. The results are shown in 
Figs. |6]-[8]for models with a flat and a flared disk, respectively. 
The behavior of the maps for the unconvolved images is as ex- 
pected and can separate the disk signal from the envelope signal 
to some extent. A ratio map of a flat disk indicates the inclina- 
tion angle. In this case the disk silhouette forms an ellipse with 
semi-major axis a (the disk radius) and semi-minor axis b, so 
the inclination is given by £ = arccos (bja). Because the radial 
extent of the modulation in the amplitude ratio map refers to 
b, the inclination can be determined. However, if the images are 
convolved with elongated PSFs, the Q- and U-maps perform dif- 
ferently, depending on the PA. of the disk. For a PA. of 0°, the 
Q-map is distorted, because of the low amplitude of the carrier in 
the region of the modulation by the disk (Fig. [8] left panel, e.g.). 
The same applies for a PA. of ~ 45° and the U-map. With a 
low-amplitude signal, PSF elongation effects become dominant. 

The maps of the amplitude ratio of the observations in J-, 
H-, and K-bands are shown in Fig. [9] In the Q-maps we find ex- 
tremal ratios in the H- and K-band at a position angle of 144° 
of the disk, so the Q-map should be less affected by elongation 
effects. The modulation in the U-maps is not significant - note 
the individual scaling of each band - and is interpreted as PSF- 
and noise effects. This is clear evidence that the disk is inclined 
and the portion in front of the star is located in the south west. 
Moreover, the radius b of the minima can be determined pre- 
cisely to 0. 12" in the H-band and in K-band (Figs.|9]and|4]right, 
dash-dotted line). If a radius of a — 0.35" and a flat disk are 
assumed, the shape of the ellipse is set and the minimum incli- 
nation angle is given by £ = arccos (b/a) s* 70°. Considering the 
height of the (flared) disk, this value is increased. This can be 
seen in Fig. [8] where the flared disk produces a greater extent as 
the flat disk for the same inclination. This results in an estimated 
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Fig. 6. Monte Carlo simulations of circumstellar envelope and disk. Left panel from left to right: Q-image, U-image, Q- and U- 
profile, Q-amplitude ratio map, and U-amplitude ratio map in H-band. The disk is dense and flat (M = O.1M , h = 0.0017?,$). The 
inclination is £ = 10°. From top to bottom: unconvolved images with P.A.= 0° for the disk, images convolved with elongated PSF, 
P.A.= 0° and convolved images with P.A.= 50°. Right panel: the same for a dense and flared disk (M = 0.1M G , h - O.OIRs). The 
profiles and maps show no significant modulation in the nearly edge-on configuration. 




Fig. 7. Same as Fig. [6] but for £ = 45°. The effects of the disk on the envelope signal is strong, as the amplitude of the profile is 
attenuated where the disk lies in front. In the maps, the PA. is clearly seen and the radial extent of the modulation refers to the 
inclination angle. 



possible inclination range of about 75°-85°, as the disk has to 
be inclined according to the signature in the ratio maps. The val- 
ues of maximum ratios increase from the J- to the K-band. This 
can be interpreted as a stronger relative attenuation of the disk 
at higher wavelength. Since the polarization is lower in K-band, 
the total attenuation might be constant. 

4.2. Modeling 

There are two main processes for producing NIR/optical po- 
larization: (i) scattering and (ii) differential absorption by non- 
spherical dust grains with the short axis primarily aligned 
along the local magnetic field, referred to as dichroic extinction 
(Bast ien & Menardlll988l) . In the case of a circumstellar enve- 
lope and disk, the first is relevant. To interpret the polarization 
data, we used a 3D radiation transfer code, based o n the Monte 
Carlo method (IWhitnev & Hartmannl Il99l 1 19931) . The model 
parameters of the presented simulations, which illustrate the sig- 
nature of the disk in the polarization images (Figs.lSfS, are de- 
scribed in the following. In addition different models for grain 
size distributions and the impacts on the polarization images are 
discussed as an attempt to explain the increasing polarization 
pattern from the J- to the K-band. 



The standard models for low-mass star formation include a 
slowly rotating, isothermal dense molecular cloud that collapses. 
The material with high angular momentum forms the circum- 
stellar disk, while the one with lower angular momentum falls 
into the proto star. We adopt ed the mode l of an infalling proto- 
stellar cloud described by iTerebev et al.l (Q~984), often referred 
to as the TSC-model. Our configuration includes both a circum- 
stellar envelope and a disk. We skipped a detailed modeling of 
bipolar cavities, sinc e we do not see any clear ev idence in the 
data. As mentioned in Whitney & Hartmannl (1 1 993b . an infalling 
envelope can easily produce the specific polarization pattern of 
reflection nebulae, which are often interpreted as scattering from 
a disk. As a consequence we expect a highly inclined disk in the 
case of 177-341 as suggested by the optical images (Pig. [5j, al- 
though we see a round, nearly centrosymmetric pattern in the 
polarimetric images. The critical parameter for reproducing the 
pattern is the mass infall rate in the TSC-model. As the enve- 
lope size is limited to about 300AU in our images (Fig. 0), the 
centrifugal radius is determined to 150AU in the model. 

The dependency of polarization intensity on the wave- 
length refers to the effective grain size in the dust envelope, 
which scatter and polari ze the emission from the central star 
(Elsasser & Stauddll978l) . The polarization spectrum and peak 
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Fig. 8. Same as Fig.|6l but fort; — 80°. Nearly seen edge-on, the disk modulation is still seen in the envelope signal. For the flat disk 
(left panel) the extent of the modulation refers to the inclination, as the height of the disk is negligible against the minor axis of the 
ellipse in projection. The extent of the highly flared disk (right panel) is dominated by the disk height. 

Table 3. Dust properties of the standard models MRN and KMH. 



Band 






MRN 








KMH 






KMHdc 




K 


ID 


g 


^rnax 


K 




g 


^max 


K 




g 


Anax 


J 


65 


0.42 


0.16 


0.81 


63 


0.46 


0.32 


0.58 


48 


0.60 


0.34 


0.58 


H 


38 


0.33 


0.06 


0.91 


38 


0.42 


0.29 


0.59 


33 


0.56 


0.31 


0.59 


K 


20 


0.21 


0.03 


0.94 


22 


0.36 


0.25 


0.60 


20 


0.50 


0.27 


0.60 



Note: The unit of k is [cnrg ']. 




Fig. 9. Plot of the test of the inclination of the disk: the ampli- 
tudes of the polarization profile (see Fig.|4| are attenuated due to 
the multiple scattering in the disk. If the disk is inclined, the 
amplitude of the part with the disk in front is more affected. 
Therefore the ratio of the amplitudes ("in front" divided by 
"behind" ) is minimal. Here the maps of the Q- (top) and U- 
amplitudes ratio in J-, H-, and K-bands (from left to right) are 
shown in an annulus of 0.1"-0.25". The position angle for the 
disk can be located around 144°. The radius of the minimum 
spots ratio along the line perpendicular to the disk can be deter- 
mined to ~0.12" in Q (H- and K-bands). 



polarization also depend strongly on the albedo oj, the opac- 
ity k, and the asymmetry parameter g of the grains. The lat- 
ter defines the relation between forward scattering and isotropic 



scattering. For the grain size distribut i on and parameters, we 
adopt ed the models f rom iMafhis et ail dl977l) . hereafter MRN 
and iKim et alj d 19941) . KMH. Table [3] gives an overview of 
the differences between the two com monly used models i n the 
NIR (data from lKim et aT] ( [19911) and I Whitney etail (fl997h ). In 
MRN, grain sizes are limited to 0.25/im, which is small com- 
pared to the wavelength of NIR-light, so the albedo and opac- 
ity drop for increasing wavelength. The KMH-model uses grain 
parameters obtained from fits to the interstellar extinction law. 
Local variations in the interstellar extinction curve are param- 
eterized with the ratio of total to selective extinction Ry = 
A(V)/E(B - V). While the standard KMH deals with the diffuse 
interstellar medium {Ry = 3.1), the alternative (here denoted as 
KM Hdc) is calculated for a dense cloud region with Ry = 5.3. 
Math is & Wallenhorstld 19811) stated that (i) an increasing R v im- 
plies an increasing size of particles, and (ii) Amax, the wave- 
length of maximum polarization, is larger for regions with Ry 
larger than the average 3.1. The peculiar extinction of 6 l Ori C 
is specified with Ry = 5.5. From the fit to the modified extinc- 
tion curve in KMH/jc, the Ry = 5.3 distribution contains signif- 
icantly fewer small- and int ermediate-sized p articles and a mod- 
est increase in larger sizes dKim et al.ll 19941) . This implies both 
a higher albedo and slightly greater g and therefore an increased 
forward scattering. The extinction is less at short wavelengths 
and flattened towards higher wavelengths. We also tested a mod- 
ified model with different grain properties only for the disks re- 
gion. If grain growth in the disk is achieved, the scattering prop- 
erties here s hould differ from the e nvelope. We adopted calcu- 
lations from iPendleton et alj d 19901) for large grains with a ~4 
times higher opacity and higher albedo, but also implemented a 
wavelength dependency, so that the scattering in K-band is in- 
creased. With high opacities and thus a high degree of multiple 
scattering, photons could be reprocessed by the disk more easily. 
Therefore more photons would be available in the outer regions 
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Table 4. Parameters used in Monte Carlo simulation. 



Parameter 


Value 


Number of photons 


10.000.000 


Stellar radius R$ 


1.5R 


Disk mass M 


0.1 M 


Outer disk radius 


150AU 


Inner disk radius 


4R S 


Inclination angle 


75° 


Position angle 


144° 


Scale height of disk h 


O.OOl^s 


Inner envelope radius 


i0R s 


Outer envelope radius 


150AU 


Mass infall rate for envelope 


1 • 10-" M Q y- 



Note: The spatial dimensions are derived from the observational data 
as described in Sect. 14. II 




[micron] 



Fig. 10. Results from simulations with different grain models in 
comparison to the observational data. The polarization extent in 
arcsec as a function of wavelength is shown (solid line, deter- 
mined as in Fig. @). The standard models MRN (dotted line) 
and KMH (dashed line) show a reversed spectral dependency. 
KMH DC (dash-dotted lines) for a dense cloud is flattened. 



of the envelope, so the scattering toward the observer could be 
enhanced. 

The parameters for the envelope and the disk are chosen as 
determined from the data in the previous Section and listed in 
Table |4] This includes the envelope and disk size, and the incli- 
nation and PA. of the disk. The mass infall rate for the envelope 
was probed in a wide range of 10 -9 — 10 6 M Q y~ l to vary the visi- 
ble extent of the polarization pattern. A value of about 10~ 9 gives 
the best results for the average radius of 0.48". Such low infall 
rates would refer to a low-mass proto star in a transition between 
late Cl ass I (embedded in e nvelopes) and Class II (with T Tauri 
disks) (IWhitnev et al.ll2003l) . 

The simulation images are convolved with the PSFs from the 
observational data to compare the results. The convolution with 
the elongated PSFs causes a change in the polarization proper- 
ties: an offset of polarization of the position angle in the direction 
of the elongation is now visible (Fig.[T3] only available electron- 
ically). This concurs with a shift in the Q- and U-profiles as seen 
in the centrosymmetric regions of the data (Fig.|4]and Sect. 14. 1 I t. 
The amount of offset depends on the specific model parameters, 
particularly the densities of disk and envelope and the disk incli- 
nation. 



4.3. Discussion 

iHennev feO'Dellld 19991) derived an inclination angle of the disk 
of 177-341 of 75°-85° by fitting photo-evapora tion model data 
to em ission line profiles of spectroscopic data. Robberto et al. 
(2002) modeled the SED for a photo-ionized proplyd and found 
that the disk emission is bright in the infrared and comparable to 
the central star. They also propose a highly inclined d isk config- 
uratio n for 177-341 by comparison to lOpm-images. Ball y et all 
(120001) interpreted HST [O III]-images of 177-341 as a disk seen 
in silhouette against the bright IF. They propose a PA. of 105°, 
although the silhouette suggests a PA. of about 135° to 145° 
because of the compact source about 0.5" east of this proplyd, 
which is likely to be an ultra compact H erbig-Haro object (177- 
341b in Fig. [TJ. IVicente & Alvesl d2005l) performed a deep anal- 
ysis on the HST data and determined the disk size from the IF 
diameter. For a distance of 0.01-0.3 pc to 6 l Ori C the FUV 
(far-ultraviolet) photons dominate the diss ociation process and 
the disk radius is typically > r IF (Johnst one et al.l 1998; 
IStorzer & Hollenbach][l999l) . In the case of 177-341 the diam- 
eter of the IF is -0.8" or 350 AU. 

Considering these constraints on the disk and envelope con- 
figuration we find good agreement of the parameters with our 
polarimetric analysis. The compliant diameter and PA. justify 
the deconvolution process and the introduced polarization am- 
plitudes ratio maps. The inclination range obtained is also con- 
sistent with previous estimates. In the special case of a polarized 
envelope and a circumstellar disk, polarimetric differential im- 
ages can be interpreted with this technique even under poor SNR 
and AO conditions. With the obtained values for the envelope 
size and the size, PA. and inclination of the disk the simulations 
show similar polarization images. 

Ac cording to a common model of the eroding disk mech- 
anism dBallvl 12003). first the outer parts of the disk are being 
destroyed. One way for disks to survive the photo-evaporation 
by UV-radiation is for the dust grains to grow to sizes with 
radii > 1cm. Evidence of grain growth in circumstellar disks in 
Orion has previously been found, because the ou ter portions of 
the gi ant disk of 1 14-426 contains grains > 2fim (Shupin g et al.l 
2003). The efficiency of grain growth is predicted to be highest 
in the center of the disk where the highest densities and tem- 
peratur es are found and p hoto-evaporation does not operate effi- 
ciently dThroop et al.ll200ll) . 

However our comparison of different dust models cannot ex- 
plain the peculiar spectral polarization of 177-341. The models 
for dense clouds with high Ry result in improved scattering at 
longer wavelengths, but do not reproduce the increasing polar- 
ization pattern extent from the J- to the K-band (Fig.fTOli. Also a 
simple model (KMHdc+) of large grains in the disk shows only a 
moderate increase in the extent. As a conclusion, advanced mod- 
eling is required including significantly different dust properties. 
In particular the opacity distribution within the envelope might 
be worth a de t ailed i nvestigation in future model improvements. 
Hughes et al. (2008) recently showed discrepancies between the 
extent of dust continuum emission and molecular gas emission, 
suggesting a different shape for the edge of the disk. An alterna- 
tive disk model with a tapered edge agrees more with observa- 
tions than the commonly used truncated power-law model. 

5. Summary and conclusions 

We have presented polarized emission maps of the circumstellar 
dust distribution of the giant proplyd 177-341 and found a com- 
plex morphology including both a circumstellar envelope and an 
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Fig. 11. Schematic overview of our model and results for the 
proplyds 177-341 and 178-342. 177-341 has an inclined disk and 
an illuminated dust envelope. The diameter, inclination, and P.A. 
of the disk and the size of the envelope are drawn according to 
the results of the polarimetric analysis. 

accretion disk. The size of the disk was estimated with the diam- 
eter of the IF and confirmed in the polarization degree images. 
We extracted the inclination of the disk from a detailed study 
of the polarization profiles in the Stokes images. For the P.A. 
of the disk, we presented evidence of part of the disk seen in 
silhouette in the optical images, and we find matching position 
angles in the polarization amplitude ratio and in the deconvolved 
polarization-degree images due to multiple scattering in the disk. 
An overview of the results of the envelope and disk configura- 
tion is given in Fig.QT| The size of the dust envelope is obtained 
from the polarimetric differential images (Fig. and increases 
from the J- (blue) to the K-band (red). In the H-band it is compa- 
rable to the diameter of the ionization fronts head in the optical 
(Fig. Q] right). 178-342 has no resolvable dust distribution, as 
the images in Fig. [3] show and modulates the shape of the IF 
of 177-341 in the optical (Fig.[TJ. The wavelength dependency 
of polarization intensity and the extent of the polarization pat- 
tern cannot be reproduced in the models by large grains in the 
circumstellar material. The comparison of common dust models 
shows evidence of the contribution of larger grain to polarization 
and scattering effects, but another mechanism must be present to 
explain the observational results. 

We showed with our 4 m-class data that the technique of 
ground-based polarimetric differential imaging is suitable for 
mapping the dust emission very close to the central star (here 
about 0.1"). This is the case for the largest and brightest pro- 
plyds. Other proplyds close to the Trapezium might be worth ob- 
serving with an 8 m-class telescope. The proplyds 158-326 and 
166-316 are the most suited candidates, as these show extended 
emission and a modulated trefoil pattern in the polarimetric dif- 
ferential images. With the introduced PDI analysis techniques, 
we are able to constrain disk parameters, such as inclination, 
PA., and size even when the disk structure itself is not resolv- 
able. 
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Fig. 12. Comparison of the different grain models used in the simulations. From left to right: MRN, KMH, KMH/x; an d KMHoc+ 
large grains (models explained in the text) in the disk region (each in J-, H-, and K-bands from top to bottom). The polarization 
vectors and the extent of the Q- and U-patterns are shown, determined as described in Sect. 14.11 The images have been convolved 
with the PSF of the observations. The deviation from circular symmetry of the vectors in the outer regions is an effect of the 
elongated PSF. 
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Fig. 13. Polarization vector maps before (top) and after shifting the Q- and U-images (bottom). The offset is determined from the 
asymmetry around zero for the azimuthal Q- and U-profiles close to the star (0.35" in J-band). The shifting compensates an offset in 
polarized flux with fixed angle in the corresponding aperture induced by the elongation of the respective field PSFs in the direction 
of 6 l Ori C. The offset dominates in the outer regions, where the polarization degree is high, but the total polarized flux is low. As 
expected, the effect is irrelevant close to the star. However, in H-band (middle) the effect is very strong for the reflection nebula and 
is interpreted to be dominated by errors in the polarimetric calibration process, as the flux in the outer regions is close to the noise 
and instrumental polarization. Hence, the shifting based on the higher flux close to the star provides a bootstrap calibration since 
the obtained vector pattern is what we expect for the reflection nebula. 



S. Rost et al.: Near-infrared polarization images of the Orion proplyds, Online Material p 4 




a? : i it.-- o.s &fi cu i>i n,4 ojs d.g 



Fig. 14. Effects of convolution with an elongated PSF on polarization. The model is simulated in H-band and includes a flat disk 
with moderate mass (M = 0.05M Q , h = 0.001/?s) and an inclination of £ = 75°. The images were convolved with artificial PSFs 
of 0.3" FWHM with varying eccentricities and rotated from 0° to 180°. Top: integrated polarization degree and angle in an circular 
aperture of 0.3". The values are given in relation to the results for a circular symmetric PSF (Pq, 6>o). Bottom: shift of the Q- and 
U-images relative to the maximum amplitude of the profile. 
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Fig. 15. PSF estimation from the FOV in comparison to the calibration star PSFs. Three point-like sources are chosen to extract the 
PSF from the field. The elongation is higher in the field PSFs, the performance in J-band is much lower. The eccentricities obtained 
by fitting the field-PSFs are 0.27, 0.35 and 0.26 for J-, H- and K-band. 



